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We present the experimental generation of large effective magnetic fields for ultracold atoms using 
photon- assisted tunneling in an optical superlattice. The underlying method does not rely on the 
internal structure of the atoms and therefore constitutes a general approach to realize widely tunable 
artificial gauge fields without the drawbacks of near-resonant optical potentials. When hopping in 
the lattice, the accumulated phase shift by an atom is equivalent to the Aharonov-Bohm phase of a 
charged particle exposed to a staggered magnetic field of large magnitude, on the order of one flux 
quantum per plaquette. We study the ground state of this system and observe that the frustration 
induced by the magnetic field can lead to a degenerate ground state for non- interacting particles. 
We provide a local measurement of the phase acquired by single particles due to photon-assisted 
tunneling. Furthermore, the quantum cyclotron orbit of single atoms in the lattice exposed to the 
effective magnetic field is directly revealed. 



I. INTRODUCTION 

The application of strong magnetic fields to two-dimen- 
sional electron gases has led to the discovery of seminal 
quantum many-body phenomena, such as the integer and 
fractional quantum Hall effect pQ. Ultracold atoms con- 
stitute a unique experimental system for studying such 
systems in a clean and well controlled environment and 
for exploring new physical regimes, not attainable in typ- 
ical condensed matter systems [2j [3]. However, charge 
neutrality of atoms prevents direct application of the 
Lorentz force with a magnetic field. An equivalent ef- 
fect can be provided by the Coriolis force in a rotating 
atomic gas, which leads to the formation of quantized 
vortices in a Bose-Einstein condensate [4 j. The regime 
of fast rotation, in which the atomic gas occupies the 
lowest Landau level, was achieved in Refs. [5] but the 
amplitude of the effective gauge field remained too small 
to enter the strongly correlated regime that requires a 
number of vortices on the order of the particle number 
[2j [6] . An alternative route consists in applying Raman 
lasers to the gas in order to realize a Berry's phase for a 
moving particle [7[ [8] . The effective gauge fields gener- 
ated in such a setup resulted in the observation of a few 
vortices, but were still far from the strong- field regime. 

Ultracold atoms in optical lattice potentials seem to 
constitute a promising system to enter the high-field 
regime. This has been shown in several theoretical pro- 
posals [SHU], which are based on the coupling between 
the motional degree of freedom and the internal structure 
of the atoms. The main challenge in lattice systems is to 
engineer complex tunneling amplitudes, that are spatially 



dependent, K = \K\e l ^ v \ This representation of the ef- 
fect of a vector potential is known as the Peierls substi- 
tution, where </>(r) is the Peierls phase [12] [13]. Recently, 
tunable Peierls phases have been realized experimentally 
in an effective ID "Zeeman lattice" using a combination 
of radio-frequency and Raman coupling [14] and in peri- 
odically driven one-dimensional lattices [T5] . 

In this Letter, we describe the creation of effective 
staggered magnetic fields for ultracold atoms in a two- 
dimensional optical lattice. The method we use is closely 
related to the proposal by Jaksch and Zoller [9] and sub- 
sequent work [10 j. However, it does not rely on the in- 
ternal structure of the atom. We make use of the well 
known phenomenon of photon-assisted tunneling [T6UT8] , 
as proposed by Kolovsky [19] , using a pair of far-detuned 
running- wave laser beams. This allows us to engineer 
spatially-dependent Peierls phases 0(r). The resulting 
effective magnetic flux through one unit cell of the lat- 
tice is fully tunable by the wavelength of the running- 
wave beams and their geometry. In our setup, the mag- 
netic flux per four-site plaquette is staggered with a zero 
mean, alternating between tt/2 and — tt/2 [20]. We study 
the nature of the ground state in this optical lattice by 
measuring its momentum distribution and show in par- 
ticular that the frustration associated with the effective 
magnetic field can lead to a degenerate ground state for 
single particles, similar to the prediction of Ref. [2T] . 
We also study the quantum cyclotron dynamics of single 
atoms restricted to a four-site plaquette. Te main results 
we present here are published in Ref. [22] . The purpose 
of this paper is to provide a more detailed discussion of 
the underlying method and to present several techniques, 
that are relevant for the experimental realization. 
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II. GENERAL IDEA 

In this section we briefly discuss the general idea of 
the method used in our experiment following mainly the 
proposal by Kolovsky [19 j. Consider atoms in a 2D lat- 
tice that is tilted and subjected to a periodic driving force 
along one direction (x-axis). In the tight-binding approx- 
imation the dynamics of this system can be described by 
the following Hamiltonian: 




H = -Jx^ ( a L+l,n«m,n + h.C.) 
m,n 

~~ (^ra,n+l^n,™ + h.C.^ 

+d x [F + F u COs(ujt + (j))] ^ m ^m,n , (1) 



where a m ^ n (a^ n ) annihilates (creates) a particle on site 
(m,n) and n m?n is the corresponding number operator. 
The strength of the tunnel coupling between neighboring 
lattice sites in the x (^-direction is given by J x ( y ) and 
d x ( y ^ is the corresponding lattice constant. The magni- 
tude of the linear potential is given by F, while the rel- 
ative driving force between neighboring sites is F^. For 
resonant driving, where uj matches the Bloch frequency 
ujb = d x F/h, and for weak coupling compared to the 
energy offset between neighboring sites 



J x < d x F, (2) 
the system can be described by an effective Hamiltonian 

m,n 

-J y J2 (^L,n+l^m,n + h.C.) , (3) 

where K = J X J{F UJ / F) is the effective tunnel coupling 
due to the driving [23[ [2^ and J{x) is the zeroth order 
Bessel function of the first kind. Experimental control 
over the phase <j) has been demonstrated in Ref. [25] . 

In order to mimic the physics of a charged particle in 
a magnetic field, the phase <\> has to depend on the po- 
sition in the lattice so that a particle that tunnels along 
the borders of one plaquette picks up a phase that can 
be interpreted as an Aharonov-Bohm phase. Experimen- 
tally this can be achieved using two running-wave beams 
with wave- vectors ki^ (|ki| ~ |ki2 1 = k) and frequencies 
co>i 5 2 such that uj\ — UO2 = ojb- The local optical potential 
created by these two running-wave laser beams is pro- 
portional to 



FIG. 1: Lattice setup to create tunable artificial magnetic 
fields. It consists of a 2D array with lattice constants d x , y 
and a linear potential Fx along the x-direction. A pair of 
running- wave beams with wave vectors |ki| c± | lc2 1 = k and 
a frequency difference cji — 002 = ojb is used to modulate the 
lattice on resonance, so that the corresponding optical poten- 
tial Vk gives rise to a spatially dependent complex coupling. 
Here a snapshot of the running- wave potential Vk with wave 
vector q is shown. The resulting effective magnetic flux <E> per 
plaquette (highlighted by the yellow box) is fully determined 
by the opening angle 20 and the frequency of the running- 
wave beams. Two examples are illustrated here for k = 7v/d y 
and (a) = 30° and (b) 6 = 15°. 



where q = (ki — k2)/2, |q| = ksiriQ and is half of the 
opening angle between the two beams. Figure [T] shows 
two examples for different geometries of the running- wave 
beams. Due to the interference between those beams the 
phase (j) increases linearly along the ^-direction according 
to Eq. Q so that on each lattice site 0(m, n) is given by 
the following relation: 



0(m, n) = 27ran, with a 



2tt jd y 



(5) 



Note that even though the photon-assisted tunnel cou- 
pling occurs along the ^-direction, the phase </>(m, n) is 
solely determined by the lattice constant d y and the pro- 
jection of the wave vector q on the ?/-axis. The value of 
the flux $ = 27ra is fully tunable through the frequency 
and the angle between the pair of running- wave beams. 
For the situation depicted in Figjlja), where 6 = 30°, 
we obtain = tt and for Fig. [TJb) , where = 15°, 

<I> w 0.527T. 

One possible scheme to realize the relative driving am- 
plitude F u along the direction of the linear potential F 
(x-axis) as described by the Hamiltonian in Eqs. ([T]) and 
(pi) would be to shift the crossing point of the running 
wave beams away from the atom position along that axis 
by an amount xq [19 . Then the local optical potential 
for the geometry depicted in Fig. [l] is given by 



■x ) cos 2 (qy - u B t/2), (6) 



Vk(t) (x cos 2 (qr — cjb^/2), 



(4) 



where V(r) is proportional to the envelope of the light in- 
tensity. This leads to a driving force F u = (dV(r)/dx)/2 
and an additional static contribution, so that the linear 
potential is given by F + (dV(r)/dx)/2. 
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III. EXPERIMENTAL SETUP 

In our experiment the 2D lattice is created by two 
standing waves of laser light at A s = 767 nm ("short 
lattice", see also Ref. [22] ). Instead of a linear potential 
in the ^-direction we superimpose another standing wave 
along that axis with twice the wavelength ( "long lattice" , 
Xi « 2A S ) to create a superlattice potential [28 , 29 , where 
every second lattice site is lifted by an energy offset A as 
depicted in Fig. [2|a). 

The lattice parameters are chosen such that A is 
much larger than the nearest-neighbor coupling J x in 
order to inhibit tunneling. By resonantly modulating 
the lattice with two running-wave beams with wave- 
length ~ 1534 nm (5uj = A/ft, with 5uj = uj\ — ^2) 
tunneling is restored with an effective coupling strength 
K = J x Vg/(2V2A), where V£ is the depth of the optical 
potential created by the pair of running- wave beams. For 
a detailed derivation of this relation see Sect. IIVBI To 
avoid cross-interference between the long lattice and the 
pair of running-wave laser beams we offset them in fre- 
quency by 160 MHz. As explained in the previous sec- 
tion, the modulation by the running- wave potential leads 
to a complex phase factor that depends on the position 
in the lattice and is determined by the geometry and the 
wavelength of the running- wave beams. In this setup the 
wave vectors are chosen to be |ki| c± | lc2 1 = fcs/2, with 
k s = 27r/A s , and the angle between the beams is 26 = 90° 
as illustrated in Fig. [5|b), so that the optical potential is 
given by 



V K (v) = V* cos 2 (^(e x + e y )r - - u 2 )t/2 (7) 



and q y = fc s /4. 

Let us consider an atom tunneling from a site of low 
energy (even sites) to a site of high energy (odd sites). 
Assuming that uo\ > 002, the phase factor is given by 
0(ra, n) = |(m + n) as illustrated in Fig. [2^b) and a = 
1/4. For atoms that tunnel from odd to even sites, the 
role of the two beams is reversed and the sign of the phase 
changes: 0(ra, n) = — f(m + n). The resulting spatial 
phase pattern for our system is shown in Fig. [2|c). Due 
to the staggering we realize an effective magnetic flux 
$ = 7r/2 alternating in sign along the direction of the 
superlattice potential and the system can be described 
by the effective Hamiltonian 



-J ( a m,n+l«m,n + h.C.) , (8) 



where J is the coupling strength along y. The sign of 
the phase factor is positive for even sites and negative 
for odd ones. 
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FIG. 2: Experimental setup, (a) Superposition of two stand- 
ing waves at A s = 767 nm (blue) and Xi = 1534 nm (red) lead- 
ing to a staggered energy offset. Black filled circles illustrate 
sites of low energy (even) and white circles the ones of high en- 
ergy (odd), (b) The experimental setup consists of a 2D array 
of ID potential tubes with spacing | | = \d y \ = d = X s /2. 
While bare tunneling occurs along the y direction with am- 
plitude J, it is inhibited along x due to the superlattice po- 
tential. A pair of far-detuned running-wave beams creates 
an additional running- wave potential Vk with wave vector q, 
which induces a complex tunnel coupling of magnitude K, 
whose phase depends on position. This realizes an effective 
flux ±$ per plaquette with alternating sign along x. (c) Spa- 
tial distribution of the complex phase of the effective tunnel 
coupling realized in the experiment. The shaded area high- 
lights the magnetic unit cell. Figure adapted from Ref. [22 . 



The magnetic unit cell of a staggered flux lattice con- 
sists of two plaquettes with opposite sign of the flux 
as illustrated by the shaded area in Fig. |2|c), see also 
Ref. [21 . It contains two non-equivalent sites (even 
and odd), therefore the band splits into two sub-bands 
[26j [27]. Considering the magnetic translation symme- 
tries we can make the following ansatz for the wave func- 
tion using the gauge realized in our setup: 



^ ^ ^i{m-k x d-\-n-k y d) 



ijj e m even 

ifj ( m + n ) m odd 



(9) 

with the quasi- momentum k in the ranges —7r/(2d) < 
k x < 7r/(2d) and —ir/d < k y < n/d. By inserting this 
ansatz into the Schrodinger equation using the Hamilto- 
nian given in Eq. ([8| one obtains a simple two-dimension- 
al eigenvalue equation for the amplitudes ^e(o) on even 
(odd) lattice sites: 



(10) 



with 



rr _ f — 2J cos(kyd) -K(e ik * d - ie~ ik * d ) \ 
~ V -K(e~ ik * d + ie i/c * d ) 2 J sin(fe 3/ d) y ' 

(11) 

The solutions of this equation determine the dispersion 
relation of the system for any ratio J/K. Due to the 
intrinsic structure of the wave function given in Eq.([9| 
every momentum eigenstate k has two momentum com- 
ponents at (k x ,k y ) (even sites) and (k x ,k y ) + (k s /2,k s /2) 
(odd sites). This means that for a non-degenerate ground 
state we expect to see two momentum peaks within the 
first magnetic Brillouin zone and replicas at momen- 
tum components shifted by the reciprocal lattice vectors 
(fc s ,0) and (0,2fc s ). 

IV. CALIBRATION OF THE SETUP 

In this section we give details about the underlying 
setup relevant for the experimental realization of the 
staggered flux lattice and calibration measurements to 
study the modulation induced coupling strength K. 

A. Staggered energy offset 

It is important for the implementation of the staggered 
flux lattice that the coupling between even and odd sites 
is equal and does not depend on the lattice site. An 
example of such a potential is shown in Fig. j4^a) . The 
particular shape of the superlattice structure depends on 
the relative phase ip between the short- and long-lattice 
standing waves and their depths, which can be controlled 
independently. The combined potential is given by the 
following equation: 

V{x) = V? sm 2 (k lX + ip/2) + V s x sin 2 (2k l x + tt/2), (12) 

where ki = 2tt/\i and is the depth of the long 

(short) lattice. The phase of both standing waves is 
fixed at the position of the retro-reflecting mirror. Due 
to a small frequency difference between the two stand- 
ing waves 5u = 2v\ — v s , where vi( s ) is the frequency of 
the long (short) lattice, they accumulate a relative phase 
(f oc SuL that depends on the distance L between the 
retro-reflecting mirror and the atoms, see Fig.[3^a). Over 
the extent of the atom cloud the phase can be assumed 
constant to very good approximation. For our experi- 
mental setup it changes by < lmrad over the cloud size. 
By changing the frequency of one of the beams we have 
full control over the superlattice potential. 

In order to calibrate the relative phase (/?, we mea- 
sured the energy offset A between neighboring sites using 
photon-assisted tunneling. The experimental sequence 
started by adiabatically loading a condensate of 87 Rb 
atoms without discernible thermal fraction into the 2D 
lattice for a given value of (/?, such that initially all atoms 
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FIG. 3: Calibration of the superlattice potential, (a) Illus- 
tration of the phase accumulation between the two stand- 
ing waves, long (red) and short (blue), until they reach the 
atom position. Note that the distances shown here are not 
to scale, (b) Fraction of atoms transferred to odd lattice sites 
riodd as a function of the modulation frequency 5uj — uj\ — UJ2 
for (p — 0.397T. The solid line shows the fit of a Lorentzian 
fuction to our data, (c) Energy offset A between even and 
odd lattice sites measured as a function of the relative phase 
if using photon-assited tunneling. The maximum offset cor- 
responds to (p = 7r/2, which is the staggered configuration. 
The solid line is a guide to the eye. 



populated the wells of lower energy. The lattice param- 
eters were V? = 5E rU V s x = 9 E rs and V s y = 20E rs , 
where E ri = ft 2 /(2mAf), i = Z,s, is the respective re- 
coil energy. In addition we used a third lattice along the 
vertical direction, with X z = 844 nm and V z = 20E rz , 
where E rz = h 2 /(2m\ 2 z ), to avoid coupling to longitudi- 
nal modes. Then we switched on the running- wave beams 
for a finite time, varying their frequency difference 5uj 
for different experimental runs. When Suj = A/ft we ob- 
served resonant particle transfer between even and odd 
lattice sites. Figure [3|c) shows the results of the offset 
calibration as a function of the relative phase. The high- 
est energy difference occurs for cp = tt/2 and corresponds 
to the staggered potential, where the value of the offset 
is A/h = 4.4(1) kHz. An exemplary resonance curve is 
shown in Fig. [3Fb), where we measured the fraction of 
atoms on odd sites n Q dd = ^odd/C^even + ^odd) as a 
function of the modulation frequency, here iV odd ( even ) is 
the atom number on odd (even) lattice sites. The mod- 
ulation induced coupling strength K/h for the measure- 
ments presented in this section is on the order of 30 Hz, 
while the damping time of the tunnel oscillations due to 
inhomogeneities in our system is typically on the order of 
a few ms. Therefore the amplitude of the resonance is at 
maximum 0.5, which corresponds to an equal population 
on even and odd lattice sites. To measure the atom pop- 
ulation in even and odd sites we employed a detection 
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FIG. 4: Calibration of the modulated tunnel coupling, (a) 
Schematics of a double well potential with energy offset A, 
inner- well coupling J x and inter- well coupling J' x . For these 
measurements J x was about two orders of magnitude bigger 
than J' x . (b) Modulation potential Vk-(t) as a function of time 
for even (black) and odd sites (gray). The relative modulation 
amplitude SVk is given by their difference and is shown in the 
lower graph (red curve), (c) Measured coupling strength K /h 
as a function of the inner- well coupling J x . The black curve is 
a linear fit to the data. The inset shows the tunnel oscillations 
between even and odd sites as a function of the modulation 
time tmod for J x /h — 0.77(4) kHz. The black line shows the 
fit of a damped sine- wave and gives K/h = 0.23(1) kHz. The 
horizontal error bars show the uncertainty in the calibration 
of the lattice depths and the vertical ones show the errors 
of the fit of a damped sine- wave. A typical image after the 
band mapping sequence for site-resolved detection is shown in 
the lower right corner together with a schematic drawing of 
the Brillouin zones (BZ). Atoms that are located on even sites 
appear in the first BZ while atoms in odd sites are transferred 
to the third BZ. 



sequence explained in Sect. IV B| (see also Fig. [5|. 



B. Modulation induced tunnel coupling 

In order to calibrate the modulation induced coupling 
strength K, we performed measurements in isolated dou- 
ble wells, where the coupling between even and odd sites 
is much higher than the one between odd and even sites 
(J x ^> J' x ), so that all dynamics were restricted to two 
sites only (see Fig. Qa)). The corresponding parameter 
regime for the superlattice potential is V* ^> and 
< cp < 7r/4. The magnitude of the modulated coupling 
strength is determined by the inner- well coupling J X1 the 
energy offset A and the relative modulation amplitude 



5Vk (see also Sect.JTT]). Due to the form of the optical po- 
tential created by the pair of running-wave beams given 
in Eq.([7|), there is a relative phase between the modula- 
tion of even and odd sites equal to 7r/2, so that 



V K ,even(t) = V£ cos 2 (Su • t/2) and 
V KM d(t) = ^cos 2 (^/4 + ^-t/2). 



(13) 



The relative modulation amplitude 5Vk{P) is then given 
by their difference, which can be written as 



K,odd 



(t) - Vi 



Ka 



v° 



1 (t) = sin(fc • t + est), (14) 
v2 



leading to a relative modulation amplitude of SV^ = 
Vg/y/2 between neighboring sites, see also Fig. 2[b). 
This gives the following expression for the effective cou- 
pling [H [24] : 



K = J X J 



V2 a) 



(15) 



and in the limit of A > V£, K w J x Vg/(2V2A). In 
order to verify this relation we performed a series of mea- 
surements to study the magnitude of K as a function of 
J x for A and constant. 

The experimental sequence started by adiabatically 
loading all atoms in a 3D lattice of depths V t x = 35 E r \ , 
V* = 30E rs and V z = 30E rz within 200 ms. Then we 
loaded all atoms adiabatically on the lower energy sites 
by ramping up the short lattice along x in 10 ms. After 
switching on the pair of running-wave beams for a vari- 
able amount of time tmod we measured the fraction of 
atoms in odd sites n dd> see inset of Fig.[4^c). The evolu- 
tion starts at n dd = 0, because all atoms occupy the even 
sites of the system, then they start undergoing tunnel os- 
cillations between even and odd sites, whose frequency is 
equal to 2if/h. The strength of the inner- well coupling 
J x was controlled by the short lattice depth V£ , which 
was varied between 8.5 and 13.5 E rs . The correspond- 
ing values of the tunnel coupling J x were obtained from 
a band structure calculation. The energy offset A was 
chosen to be equal to 4.4 kHz. As it is shown in Fig. 4^c) 
we find a linear dependence of the coupling strength on 
the inner- well coupling with a slope of K/J x = 0.31(1). 
From an independent calibration of the potential depth 
created by the two running- wave beams V£ = 4.2(1) E r i, 
we obtain a slope V^/(2v / 2A) = 0.34(2) and taking into 
account non-linear corrections due to the Bessel func- 
tion we obtain a slope of J(Vg/ (V^A)) = 0.32(1). For 
J x /h = 0.20(1) kHz the damping time of the tunnel os- 
cillations due to inhomogeneities in our system is on the 
order of the oscillation period. This explains the larger 
uncertainty for the measured modulation induced cou- 
pling strength K/h and most likely explains the larger 
deviation from the linear curve, see Fig. [4^c). 
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FIG. 5: Site-resolved detection sequence. Atoms initially lo- 
cated in odd sites (i) are transferred non-adiabatically to the 
third energy level of the system (ii). After ramping down 
the short lattice adiabatically the occupation of the different 
Bloch bands is preserved (iii) and can be measured using a 
band-mapping technique. 



To detect the population on different sites we first 
ramped up the short lattice to freeze all tunneling dy- 
namics and preserve the atom population on each lattice 
site. Then we transferred atoms that are located on odd 
sites non-adiabatically to the third energy level in the sys- 
tem by increasing the energy offset between neighboring 
sites. Ramping down the short lattice adiabatically pre- 
served the population in the different Bloch bands, which 
were detected by a following band-mapping sequence (see 
also Fig. [5]) [28l|29]. In the lower right corner of Fig.[4^c) 
we show a typical image obtained after the detection se- 
quence together with the Brillouin zones. The size of the 
zones is reduced along the x-direction because the short 
lattice is ramped down earlier so that we are left with 
the long lattice only before the band-mapping sequence, 
which has a lattice constant that is twice the one of the 
short lattice. 



V. GROUND STATE FOR ISOTROPIC 
COUPLING J — K 

The first measurement we present in this section was 
performed to study the ground state of the staggered 
flux lattice from its momentum distribution for isotropic 
coupling strengths J = K [22 . 

The experimental sequence started by loading a con- 
densate of about 5 x 10 4 atoms with no discernible ther- 
mal fraction into a 2D lattice within 160 ms, forming 
an array of ID Bose gases. The depth of the lattice 
was adjusted to Vf = = 14E rs , which corresponds 
to an equal tunnel coupling along both directions of 
J x /h = J/h = 31(2) Hz. Then we ramped up the long- 
lattice in 0.7 ms to V t x = 5 E r \ with tp = tt/2 and subse- 
quently decreased the short lattice in the x-direction to 
9 E rs within 1 ms. The staggered energy offset between 
neighboring sites was A/h = 4.4(1) kHz and was cali- 
brated using the measurements presented in Sect. |IV A 



The bare tunnel coupling along x was J x /h = 94(4) Hz 
and therefore much smaller than A. This results in 
an inhibition of spontaneous tunneling along this di- 
rection, which was then restored by switching on the 
running- wave beams on resonance, 8uj = 27rx4.4kHz and 
V£ = 4.2(1) E rl . This yields a value of K/h = 30(2) Hz 
and was obtained using the calibration measurements 
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FIG. 6: Ground state of the staggered flux lattice for isotropic 
coupling J — K. (a) - (b) Momentum distribution mea- 
sured after a time of flight of £tof = 20 ms for (a) simple cu- 
bic lattice without modulated coupling by the running-w ave 
beams and (b) staggered flux lattice as explained in Sect. Ill 
for J/K — 1.0(1). The theoretical momentum distributions 
are obtained by an exact diagonalization of the Hamiltonian 
m Eq. (H) on a 31 x 31 lattice with harmonic confinement. 
The red square indicates the magnetic Brillouin zone and the 
cross marks the center, (c) Dispersion relation of the lowest 
Bloch band, calculated in the tight -bin ding approximation 
using the Hamiltonian given in Eq. (11) plotted for the first 
magnetic Brillouin zone. The red crosses mark the position 
of the ground-state momentum eigenstate together with the 
one shifted due to the intrinsic structure of the wave function 
given in Eq.Q. (d) Spatial distribution of the phase for the 
ground-state wave function. The vortices of different chiral- 
ity are illustrated by the rotation of the white arrows. Figure 
adapted from Ref. [22] . 



presented in Sect. IV B After holding the atoms in this 
configuration for 10 ms, we observed a restored phase co- 
herence as shown in Fig. ^b). Most likely, this can be 
attributed to a redistribution of entropy present in the 
random phases between the ID condensates along their 
longitudinal modes. 

In Fig. [6]^a) we show the momentum distribution of 
a simple cubic lattice with lattice constant d = A s /2, 
without running- wave beams. In contrast to this sim- 
ple structure with a single quasi-momentum component 
and correspondingly single momentum peaks separated 
by 2Hk s , we observe four momentum components within 
the first Brillouin zone of the short lattice in the presence 
of the staggered flux, see Fig. [6jb). In order to under- 
stand that distribution we calculated the dispersion re- 
lation using Eqs. ( [To| ) and (11) for J = K, which is plot- 
ted in Fig. [6^c). It shows a single minimum located at 
k = (— fc s /4, — k s /4) and the eigenvector has equal weight 
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on even and odd sites, |^ e | = \i/j \. Therefore we expect 
to see two momentum components within the first mag- 
netic Brillouin zone, one at (— fc s /4, — k s /4) (even sites) 
and a second one shifted due to the additional phase fac- 
tor in the wave function for odd sites at (fc s /4, fc s /4). 
The magnetic unit cell, illustrated by the shaded area 
in Fig. |2|c), contains two lattice sites, therefore the size 
of the magnetic Brillouin zone is reduced along the re- 
direct ion compared to the one of the square lattice. All 
other momentum peaks are translated by the reciprocal 
lattice vectors (fc s ,0) and (0, 2fc s ). The measured posi- 
tions are in good agreement with the quasimomenta of 
the Bloch states of lowest energy. 

The frustration introduced by the position-dependent 
phase factors </>(m, n) causes the phase of the atomic wave 
function to be non- uniform, see Fig. [6^d). The direction 
of the arrows on each lattice site indicates the phase of 
the wave function, while their length corresponds to the 
magnitude of the local atomic density. To make the ap- 
pearance of vortices more clear we extrapolate the value 
of the phase between lattice sites and encode it in the 
background color of this plot. As one expects for an ef- 
fective magnetic flux of = 7r/2, which is one quarter 
of a flux quantum, there is one vortex per four plaque- 
ttes along the ^-direction and due to the staggering, the 
chirality of the vortices alternates along x. The atomic 
density of the ground-state wave function is uniform, il- 
lustrated by the length of the arrows in Fig. |6jd), as ex- 
pected from the solution of the eigenvalue equation, that 
predicts equal weights for even and odd sites \ip e \ — \^) \. 



VI. GROUND STATE FOR ANISOTROPIC 
COUPLING 

To further study the nature of the ground state in our 
system we varied the ratio J/K and observed the result- 
ing momentum distribution of the ground state [22] . The 
experimental sequence was similar to the one described in 
the previous section. In order to change the coupling ra- 
tio we kept the modulated coupling strength K constant 
and changed the lattice depth along y to = 10E rs , 
which leads to J/h = 75(3) Hz and J/K = 2.5(2). 

As can be seen in Fig. (7^a), we observed two addi- 
tional peaks in the momentum distribution within the 
first magnetic Brillouin zone for these parameters. This 
is due to a two-fold degenerate ground state, in agree- 
ment with the solutions of the eigenvalue equation ( flOj ) . 
The corresponding dispersion relation is plotted in Fig- 
ure[7^c). It shows two minima at momenta Icl — (—0.25 
fc s ,-0.45 k 8 ) and k n ~ (-0.25 k a , -0.052 k a ), split a- 
round k y = — k s /4 by Ak y ~ 0.40 k s . The relative 
atomic density on even and odd sites is given by the 
corresponding eigenvectors. For the eigenstate kL we ob- 
tain a relative weight |^ | 2 /|^ e | 2 ~ 6.1, which means that 
the momentum component shifted to kL + (k 3 /2, k s /2) ~ 
(0.25 k s , 0.052 k s ) (odd sites) is about six times more in- 
tense than the one at kL — (— 0.25 fc s , —0.45 k s ) (even 




FIG. 7: Ground state of the staggered flux lattice for J = 
2.5K. (a) Momentum distribution measured after £tof = 
20 ms for J/K — 2.5(2). The theoretical momentum distribu- 
tions are obtained by an exact diagonalization of the Hamilto- 
nian in Eq. |8| ona31x31 lattice with harmonic confinement. 
The red square indicates the magnetic Brillouin zone and the 
cross marks the center, (b) Histogram of the measured frac- 
tion of atoms in peaks corresponding to the lower momentum 
state, kL and kL + (k s /2,k s /2), obtained from 172 identical 
experimental runs, (c) Dispersion relation of the lowest Bloch 
band, calculated using the Hamiltonian in Eq. (11) plotted 
for the first magnetic Brillouin zone. The red crosses mark 
the positions of the degenerate momentum ground state to- 
gether with the ones shifted due to the intrinsic structure of 
the corresponding eigenstate. (d) Spatial distribution of the 
ground- state phase distribution. The atomic density shows a 
charge density wave illustrated by the length of the arrows. 
A similar structure is obtained for the second degenerate mo- 
mentum state but shifted by one lattice site. Figure adapted 
from Ref. [22]. 



sites). For the second momentum eigenstate ku the rel- 
ative weight is exactly opposite | ?/^ 1 2 / 1 ^ e 1 2 — 1/6.1 and 
the component at ku + (k s /2,k s /2) (0.25 fc s , 0.45 k s ) 
(odd sites) is about six times less intense than the one 
at ku — (— 0.25 k S: —0.052 k s ). This is in good agree- 
ment with our experimental observations, as shown in 
Fig. 7^a). The atomic density in the lattice is illustrated 
by the color brightness and the length of the arrows in 
Fig. 7^d) showing a charge density wave. The phase of 
the atomic wave function tends to align along the direc- 
tion of the stronger coupling J, thereby frustrating the 
phase relation imposed by the Hamiltonian. As a con- 
sequence, the density in every second stripe along y is 
suppressed. Figure (7^d) shows the result for the eigen- 
state kL, where the odd sites are more populated. The 
result for ku is similar but shows a higher occupation for 
even sites. 

Contrary to the case of a triangular lattice with frus- 
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FIG. 8: Projection of the momentum eigenstates (odd sites) 
on the ^/-direction as a function of the coupling ratio J/K. 
For J/K < y/2 the peaks are located at (fc s /4, fc s /4), while 
for J/K > \[2 the peaks are split due to the emergent 
ground- state degeneracy, see insets. The solid lines corre- 
spond to the minima of the lowest Bloch band calculated in 
the tight-binding approximation using the Hamiltonian given 
in Eq. (11). Figure adapted from Ref. [22] . 



trated hopping, studied in Ref. [30 , the fraction of atoms 
in each single-particle ground state does not fluctuate, as 
shown in Fig. [7Jb) . We observe equal population in both 
states as predicted for weakly interacting [21] or finite size 
systems. In trapped systems the degeneracy is slightly 
lifted and the two lowest energy states are symmetric 
and antisymmetric superpositions of the two momentum 
states. 

We also varied the ratio J/K gradually between 0.6 
and 2.8 by adjusting the ^-lattice depth to values be- 
tween 9.5£? rs and 16.5 E rs , while the modulated cou- 
pling strength K/h = 30(2) Hz was kept constant. The 
observed momentum distribution remained unchanged 
for J/K < 1.4, while above that value we observed 
the appearance of a two-fold degenerate ground state 
(see Fig. [8|. For this measurement we focused on the 
momentum peaks that correspond to odd lattice sites 
(k x = k s /4) and plot their projection on the ^/-direction 
as a function of the ratio J/K. From the band-structure 
calculation we expect a bifurcation at J/K = y/2 (solid 
line in Fig. [8|, which agrees well with our data. In 
the limit of J ^> K the two ground states approach 
k L ~ (0.25fc s ,0) and k n ~ (0.25 k a , 0.5 k s ) with a split- 
ting of Ak y = 0.5k s , which corresponds to a state with 
relative atomic density \*p \ 2 /\^e\ 2 — and reversed for 
the other component. The nature of the bifurcation is 
identical to the one predicted in Ref. [21], where it is in- 
duced by a variation of the magnetic flux amplitude at a 
fixed value of J/K = 1. 

To conclude this section we point out that there are 
two Dirac cones within the first Brillouin zone at mo- 
menta k m ~ (0.25 k a , -0.25 k s ) and k D2 ~ (0.25 k s , 0.75 
k s ), see Fig. [9] It is interesting to mention, that their 
positions remain unchanged for different ratios of J/K 
even in the regime, where there is a two-fold degener- 
ate ground state whose momentum components are split 
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FIG. 9: Band structure for (a) J/K = 1 and (b) J/K = 2.5. 
The positions of the Dirac cones are marked with red circles. 
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VII. LOCAL OBSERVATION OF THE 
EFFECTIVE MAGNETIC FIELD 

In order to probe the structure of the Hamiltonian in 
a more direct way, we performed a second series of ex- 
periments to study the local effect of the artificial gauge 
field on the level of a four-site square plaquette [22 . The 
experimental setup is shown in Fig. [lC^a). We use two 
superlattice potentials along the x- and ^/-direction to iso- 
late plaquettes so that all dynamics is restricted to four 
sites without any coupling between plaquettes, similar 
to the setup described in Sect. |IVB[ where all dynamics 
was restricted to two sites only. This allows us to isolate 
plaquettes with equal sign of the flux. The relative phase 
between the two standing waves along y was chosen to 
be (p y = 0, which corresponds to a symmetric double- well 
configuration without energy offset between neighboring 
sites. Along the x-direction the relative phase was in the 
regime < cp x < 7r/4 in order to create a tilt between 
even and odd sites. To avoid coupling to axial modes 
along the potential tubes, an additional lattice along the 
^-direction was used. The four sites of a single plaquette 
are denoted as A, B, C and D (see Fig. [lO^b)). 

The experimental sequence started by loading a con- 
densate into a 3D lattice of depth V t x = 35E r i, V t y = 
35E r i and V z = 30E rz . After that, a filtering sequence 
was applied to have at most one atom per lattice site [29 . 
By ramping up the short-lattice along x to V s x = 5 E rs 
within 10 ms and ip = 0.20(1) rad we realized a double 
well potential with energy offset A/h = 6.0(1) kHz, where 
single atoms were localized on one side. Finally the short- 
lattice along y was ramped up to = 14 E rs within 1 ms 
to create a plaquette potential, where each atom is in the 
initial state = (\A) + \D))/y/2 (J/h = 0.17(2) kHz 
and J x /h = 2.0(1) kHz). Subsequently the running-wave 
beams were switched on with duo = A/h in order to in- 
duce resonant coupling to the B and C sites. The in- 
duced coupling K/h = 0.32(1) kHz was measured inde- 
pendently. 

The following rather complex phase evolution can be 
understood in the limit of J <C K, where the dynam- 
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FIG. 10: Local observation of the phase evolution in isolated 
plaquettes. (a) Schematic illustration of the experimental 
setup, (b) The initial state is a single delocalized atom over 
the sites A and D (gray area), that is coupled to B and C 
sites using the pair of running- wave beams. On the left we 
show a typical double-slit pattern for the initial state obtained 
after £tof = 20 ms. (c) Phase evolution of the double-slit 
pattern along y (integrated along x), as a function of time 
for hw — A (blue) and hw — — A (gray). The inset shows 
the Fourier transformation for hw — A depicting two fre- 
quency components at 0.24(6) kHz and 0.62(13) kHz, in good 
agreement with theory (vertical lines). Figure adapted from 
Ref. [22]. 



ics along y would be suppressed and the initial state 
couples to the state \ip2) = (\B) + i\C))/y/2, due 
to the relative phase imprinted by the pair of running- 
wave beams. In our case the ratio is J/K « 0.5 leading 
to a more complex evolution. We measured the value of 
the phase as a function of time from the shape of the 
momentum distribution after time- of- flight. The method 
can be explained considering the double-slit interference 
obtained for an isolated double well along the ^-direction, 
where we denote the two sites as \A) and \D), and we 
consider an arbitrary single-particle state 
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FIG. 11: (a) Schematic of the four-site plaquette. (b) Bril- 
louin zones of the 2D lattice, (c) Typical momentum dis- 
tribution obtained with the two-dimensional band mapping 
sequence after 20 ms of time-of- flight. Figure adapted from 
Ref. [22]. 



termined the phase j(t) along y. 

'c 



The result is shown 
For Suj —A/ft the 



with blue data points in Fig. 10 
role of the two running- wave beams is reversed, W\ < 002 
and the spatial phase pattern induced by the running- 
wave beams changes sign. Therefore we observed a sign 
reversal for the phase evolution. 

We also observed the dynamics in real space in order 
to exhibit the influence of the artificial gauge field on the 
particle flow. The experimental sequence was the same 
than the one described above, only the depths of the 
short lattices were changed to Vf = 7 E rs and V]* = 
10 E rs in order to obtain different coupling strengths 
J x /h = 1.2(1) kHz and J/h = 0.50(2) kHz. By gener- 
alizing the site-resolved detection sequence discussed in 
Sect. |IVB| for a system of isolated double wells to pla- 
quettes, we measured the atom population on each site 
N q (q = A, £?,C, D). First we applied the non-adiabatic 
mapping sequence along x and y in order to transfer the 
populations N q to different Bloch bands. Atom popula- 
tion originally in B sites is therefore transferred to the 
third energy level in the x-direction, while population in 
D sites is transferred to the third level along y. A subse- 
quent band-mapping technique allowed us then to deter- 
mine the occupation of different Bloch bands by count- 
ing the atom numbers in different Brillouin zones (see 
Fig. [Iljb)). The different colors show the connection 
between Brillouin zones and their corresponding lattice 
sites (see Fig. 11 'a)). A typical image obtained after 
20 ms of time-of- flight is shown in Fig. 11 'c). 

From these images we extracted the average atom po- 
sitions 



V2 



(\A) 



\D)) 



(16) 



with phase 7. The observed density distribution after 
a sufficiently large £tof is proportional to cos(k t y + 7) 
times an envelope determined by the Wannier function, 
where k t = mdj Mtof [28] [29]. For the initial state l^i), 
7 = and we expect a symmetric momentum distribu- 
tion. A typical image is shown in Fig. 10 ^b). For the 
following evolution in isolated plaquettes, we first inte- 
grated the momentum distribution along x and then de- 



/ y\ (-N A + N B + N c -N D )d 

{X) = 2N (17) 

m {-N A -N B + N c + N D )d 

with N being the total atom number. Initially all atoms 
occupied the A and D sites with equal weight and the 
evolution starts at (X) (t = 0) = -0.5 and (Y) (t = 
0) = 0. After switching on the pair of running- wave 
beams we observed a coherent particle flow inside the 
plaquettes towards the B and C sites (see Fig. 12 [b)). 
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FIG. 12: (a) Cyclotron orbits of the average particle posi- 
tion obtained from the mean atom positions (X) /d (b) and 
(Y) /d (c) for J/h = 0.50(2) kHz, K/h = 0.28(1) kHz and 
A/h = 4.9(1) kHz. Each data point is an average over three 
measurements. The inset in (a) shows the theoretical curve 
calculated for <j> — 0.80 x tt/2 and a 1/e— damping time of 
13 ms obtained from damped sine fits to (X) /d and (Y) /d. 
Figure adapted from Ref. [22] . 



Thereafter, the particle current showed deviations along 
the ^-direction (see Fig. 
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c)), which is reminiscent of 
the Lorentz force acting on a charged particle in a mag- 
netic field. As shown in Fig. [l2|a), the mean atom po- 
sition followed an orbit that is a small-scale quantum 
analog of the classical cyclotron orbits for charged par- 
ticles. This coherent evolution is damped due to spa- 
tial inhomogeneities in the atomic sample. Having in- 
dependently calibrated the values of J and K, we fit 
from the measured atom dynamics the value of the mag- 
netic flux (j> = 0.73(5) x 7r/2. The difference from the 
value (j) = 7r/2 expected for a homogeneous lattice stems 
from the smaller distance between lattice sites inside the 
plaquettes when separated with an additional superlat- 



tice along the ^-direction. For the parameters used in 
Fig. 12 ^a)-(c) we calculate a distance d y = 0.78(1) x A s /2 
along y yielding <ft = 0.80(1) x 7r/2, which qualitatively 
explains the measured flux value. Residual deviations 
might also be due to an angle mismatch between the two 
running- wave beams and the lattices beams. 



VIII. CONCLUSION 

In conclusion, we have demonstrated a new type of op- 
tical lattice that realizes strong effective magnetic using 
a method that does not rely on the internal structure 
of the atom and is therefore applicable to a large vari- 
ety of different systems. We have shown that the atomic 
sample relaxes to the minima of the magnetic bandstruc- 
ture, realizing an analogue of a frustrated classical spin 
system. However, the spatial average of the magnetic 
flux is zero, hence the Bloch band is topologically trivial 
[TO} [27l I5T] . By using a superlattice potential with more 
than two non- equivalent sites [10 or a linear tilt poten- 
tial [9], it is possible to create a lattice with a uniform 
and non-zero magnetic flux. This system would realize 
the Harper Hamiltonian [32] and lead to the fractal band 
structure of the Hofstadter butterfly [13]. In particular 
the lowest band would exhibit a Chern number of one and 
be analogous to the lowest Landau level [8j|9j[TTJ|33]. Our 
work constitutes an important step towards the study of 
quantum Hall effect with ultracold atomic gases and the 
creation of strongly-interacting liquids such as the Laugh- 
lin state [34 . In addition the staggered flux lattice can 
be used to measure the 7r-flux associated with a Dirac 
cone using the concept of Zak phases following a recent 
proposal by D. A. Abanin et al. [35]. This method was 
already demonstrated in the experiment for measuring 
topological invariants in ID systems [36] . 
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